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INFLUENCE  OF  MAGNETIC  SHEAR  ON  THE  COLLISIONAL  CURRENT  DRIVEN 

ION  CYCLOTRON  INSTABILITY 


I.  INTRODUCTION 

Intense  field  aligned  currents  observed  In  the  auroral  E  and  F  regions 

(Park  and  Cloutier,  1971,  Cahill  et  al.,  1974,  Kelley  et  al.,  1975)  were 

conjectured  to  be  responsible  for  the  generation  of  electrostatic  Ion 

cyclotron  waves.  Since  the  Ionosphere  Is  colllsional  in  the  E  and  F  regions 

the  field  aligned  currents  give  rise  to  negative  energy  waves  which  grow  due 

to  the  dissipation  associated  with  electron-neutral  collisions  (Chaturvedl, 

1976).  These  are  electrostatic  waves  and  lead  to  density  irregularities  in 

the  medium.  In  fact,  radar  measurements  in  the  E  region  detect 

irregularities  collocated  with  auroral  electrojets  (Greenwald  et  al.,  1975), 

and  a.c.  field  measurements  (Kelley  et  al.,  1975)  have  been  identified  as 

being  due  to  unstable  ion  cyclotron  waves  (Drummond  and  Rosenbluth,  1962  - 

collisionless  domain;  Kindel  and  Kennel,  1971  -  weakly  colllsional 

domain).  Sounding  rockets  launched  from  Syowa  station,  detected 

transversely  propagating  electrostatic  plasma  instabilities  of  scale  sizes 

around  200  km  in  association  with  strong  field  aligned  currents  (Ogawa  et 

al.,  1981).  Chaturvedl  (1976),  using  a  local  fluid  analysis,  showed  that, 

in  a  partially  ionized  plasma,  field  aligned  currents  can  support  almost 

transversely  propagating  ion  cyclotron  waves  owing  to  resistivity 

experienced  by  electron  parallel  motion  (electron-neutral  collisions),  while 

the  ion-neutral  collisions  were  found  to  have  a  mild  stabilizing  effect. 

The  field  aligned  currents  observed  are  usually  -  1-10  pA/m  (Kelley  et  al. , 

1975),  although  recently  Burke  et  al.  (1983)  reported  observations  of  large 
2 

currents  '■  135  pA/ra  .  These  field  aligned  currents  can  generate  magnetic 
shear  which  in  turn  affect  the  mode  structure.  This  is  especially  so  in  the 
case  of  collisionless  current  driven  ion  cyclotron  waves  (Ganguli  and 
Bakshi,  1982). 

In  this  brief  report  we  investigate  the  effects  of  a  self-consistent 

magnetic  shear  on  the  colllsional  current  driven  ion  cyclotron  instability 

(CIC).  We  find  that  the  magnetic  shear  has  no  noticable  stabilizing  effects 

on  the  current  driven  ion  cyclotron  instability  in  the  domains  of  interest 

pertaining  to  auroral  conditions  (characterized  by  V,  -  0.5  -  5  km/s, 

2  Q  4-1 

corresponding  to  currents  of  order  10  -  100  pA/m  ,  and  v  ^  10  s  ,  where 
Manuscript  approved  March  20, 1984. 


is  the  electron  drift  velocity  along  the  magnetic  field  and  is  the 

electron-neutral  collision  frequency).  This  is  because  the  self-consistent 

magnetic  shear  produced  by  the  currents  observed  in  the  auroral  region  is 

very  small.  This  result  is  in  contrast  to  the  major  effects  found  by 

Ganguli  and  Bakshi  (1982)  in  the  collisionless  domain.  Their  kinetic 

treatment  shows  that,  due  to  the  nonlocal  boundary  conditions,  magnetic 

shear  affects  the  collisionless  current  driven  ion  cyclotron  instability  in 

two  ways:  One,  due  to  an  explicit  shear  dependent  term  of  0(p  /L  )  which 

s  3 

vanishes  for  p  /L  -*•  0.  Two,  due  to  a  shear  independent  term  containing  the 
s  s 

derivatives  of  the  potential  function  defined  in  equation  18.  In  the  fluid 

example,  the  second  term  happens  to  be  zero  (since  the  derivative  of  Q, 

defined  in  Eq.  18,  vanishes  at  Xq  which  is  defined  in  Eq.  19)  and  hence  the 

lack  of  the  significant  effects  in  the  p  /L  ♦  0  limit.  Furthermore, 

s  s 

for  p  /L  «  1  magnetic  shear  localizes  the  wave  packet  around  k  /k  «  1, 
s  s  z  y 

indicating  that  the  mode  is  almost  perpendicularly  propagating,  and  stronger 
magnetic  shear  stabilizes  the  instability.  The  paper  is  organized  as 
follows.  In  the  next  section  we  give  the  theory  and  derive  the  nonlocal 
mode  structure  equation.  In  the  third  section  we  present  the  results,  and  in 
the  last  section  we  apply  the  results  to  the  auroral  ionosphere  and  discuss 
future  work. 

II.  THEORY 

We  consider  a  partially  ionized  plasma  such  as  the  one  encountered  in 
the  auroral  E  region.  The  geometry  used  in  this  analysis  is  as  follows. 
The  magnetic  field  is  aligned  with  the  z-direction  along  which  flow  the 
currents  that  drive  the  CIC.  Because  of  the  self-consistent  shear  produced 
by  the  field  aligned  currents,  the  magnetic  field  acquires  a  component  along 
the  y-direction  as  x  is  varied, 

a  a 

B(x)  *  B  (£  +  x/Lg  2)  •  ( 1 ) 

This  implies  that  the  parallel  wave  vector  becomes  a  function  of  x, 

k  =*  k°  +  k  ( x/L  j  ,  (2) 

z  z  y''  s' 


2 


where  Lg  is  the  shear  length.  The  self-consistent  shear  is  calculated  from 
the  Maxwell's  equation 


7  x  i  -  r  ± "  f- 


(3) 


where  is  the  drift  velocity  parallel  to  the  z-direction,  and  is  given  by 


-1 


Ls  ■  0  P.le./V 


(4) 


1/2 

where  cg  *  (T^/m^)  ,  pg  ■  cg/fl^,  is  the  ion  cyclotron  frequency,  and 


a  =  Us/Ci2  <M/mHu)pe/ne-* 


(5) 


where  o>  and  fl  are  the  electron  plasma  and  cyclotron  frequencies, 

P  ®  i 

respectively,  c  is  the  speed  of  light,  and  M  and  m  are  the  ion  (NO  )  and 

electron  masses,  respectively. 

The  basic  equations  describing  the  problem  are  as  follows.  The 
electron  and  ion  continuity  equations  are  given  by 


3n 


-  0 


3t  v~a-a' 

where  a  indicates  the  species.  The  equation  of  motion  for  the  ions  is. 


(6) 


V.  *  B 


Mn 


^  ^  -  VPi  +  en(-V*  +  — - J  -  Mnvi  , 


(7) 


while  the  equation  of  motion  for  electrons  is. 


V  x  B 

0  3  -  Vp  -  en(-V$  +  — - •)  -  mnv  V  , 

e  v  c  J  6 


(8) 


where  v  and  v  are  electron-neutral  and  ion-neutral  collision  frequencies, 
e  l 

respectively,  pa  »  na^a  *s  c^e  pressure  of  species  a,  and  d/dt  * 

3/3t  +  V/V.  In  this  paper  we  ignore  the  electron  inertia,  and  assume  that 
ve/&e  «  1*  and  «  I.  We  also  Ignore  the  viscosity  and  thermal 

conductivity  terms  because  these  terms  only  introduce  a  numerical  constant 
that  multiplies  the  growth  rate  and  do  not  introduce  any  additional  k2/ky 
terms  (see,  for  example,  Chaturvedi  and  Kaw,  1975).  We  also  ignore  the 
elctromagnetic  effects  and  consider  only  to  electrostatic  perturbations.  We 
derive  the  mode  structure  equation  assuming  that  the  perturbed  quantities 

3 


vary  as  f  '  f(x)  exp  -  k  y  -  k  zJJ  .  We  note  that  we  do  not  Fourier 

y  z 

analyze  in  the  x  direction  since  the  magnetic  field  is  assumed  to  vary  in 
that  direction. 

From  Eq.  (8),  we  obtain  the  perpendicular  and  parallel  (to  B) 
components  of  the  electron  velocity 


V  -i- 

-■  n 


e  (,l+v|/Q|J 
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7.P 

1  e  e 


mn 


-  5  V)  «  i 


7,p 

-  I  — t-S-  -  S-  7 


mn 


m  1T'  fl 


o1! 


i  7  P 

V  —  [ — —  -  —  V  <pj  . 

ez  v  1  mn  m  z  J 


(9) 


(10) 


From  Eq.  (10)  we  define 


V  =  Vu  -  -S.  (£  V  4°) 
H  o*  Vn  w  J 


B  zT 


(11) 


From  Eq.  (7)  we  obtain  for  the  ions 

iu>/fl 


2ti" 


7,p, 


L  i+ci«-+v1J  /a  Jj 


1 1 2  /n  2 1  V  +  Q  .  Mn6-1  X  il 


and 


(12) 


iZ  -id) 
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7  P 

z*e  e 


Mn 


-SVJ 


(13) 


where  we  have  assumed  Tg  T^,  and  $-$+<)>  ($  being  the  perturbed 
potential).  Substituting  Eqs.  (9)  thru  (13)  in  Eq.  (6)  we  obtain,  after 
subtracting  the  electron  continuity  equation  from  the  ion  continuity 
equation,  the  nonlocal  mode  structure  equation 


.2,  2  _2  .  2  2 


32+  +  -Q1-kycsj  +  l(>m/MjkyVe(,(j-kzVdj  +  io^  k; 


3x 


klcl  '  iU/Mjve(u.-kzVdJ 


♦  -  0, 


(14) 


where  k  *  k  +  k  x/L  and  where  quasi-neutrality  has  been  assumed.  The 
z  z  y  s 

local  limit  of  Eq.  (14)  (Lg  -*■  •  and  3/3x  «  ky)  yields  Q(x)  -  0,  where  Q(x) 


» 


I 


I 


| 


» 


4 


l 
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are  Che  second  two  terms  in  the  above  equation.  This  gives  the  results 
obtained  previously  by  Chaturvedi  (1976).  The  linear  local  growth  rate  is 

v2  k  •  V  v 

(1S> 

2M  k  r 

z 

where 


2  2  2  2 

ur  =  flj  +  k  c  . 
r  i  y  s 


0)  =  0)  +  iY 

r 

Gq.  (15)  leads  to  the  threshold  condition 

k2  v 

t-  >«ru  *M-). 

k  v 
y  e 

Note  that  the  nonlocal  treatment  allows  for  the  formation  of  a  wave  packet 
in  the  x-direction  compared  with  the  plane  wave  [exp(ik^x)]  solutions  of 
Chaturvedi  (1976). 

We  solve  Eq.  (14)  analytically  and  numerically  to  determine  the 
eigenvalues  and  eigenfunctions.  We  follow  the  techniques  developed  by 
Ganguli  and  Bakshi  (1982)  to  solve  Eq.  (14)  as 

-4-  +  Q(xH  **  0,  (17) 

3x 


where 


[(.<u2-nj-k2c2)  +  iuv  !k2x2/L2  +  i(m/M)  v  (u-k  V,x/L  jk2 
1  v  i  y  s' _ iJ  y  s  ev  yd  s'  y 

k2c2x2/L2  -  i(ra/Ml  v  (u>-k  V,x/L  ) 
ys  s  v  *  yd  s' 


We  expand  Q(x)  around  £  =  x-x  ■  0  to  0(S  ;  where  x.  is  the  value  of  x  for 

o  ° 

which  the  local  growth  rate  is  a  maximum, 


_o_  2  m 

P  *  k  V. 
s  yd 


. • 


Now,  we  have 


+  (a  +  n2J  Kn)  -  0 

3ri 


(20) 


where  n  ■  C  +  Qq/Qq*  and  a  »  Qq  -  Q'/2Q^'’.  The  primes  indicate  the 
derivatives  with  respect  to  x  evaluated  at  x  *  x  . 

Equation  (20),  which  is  in  the  form  of  a  Weber  equation,  has  the 
solution  (for  the  l  *  0  mode) 

<Kx)  *  $  exp[-  (x  -  xmoj2/x„]  exP  Uq(x)J  (21) 


where  l  is  the  order  of  the  mode  and  q(x)  is  a  real  function.  The 
wavepacket  peaks  at  (Ganguli  and  Bakshi,  1982) 

r4q:  /Qri-Qr/2j1/2j 


mo 


x  - 
o 


Re  C~Q'"/2J 


1/2 


(22) 


The  dispersion  relation  is  given  by 


«*„>  -  <2*  +  1)  l-%-J1/2  »  /L 

o 


(23) 


It  can  be  easily  shown  from  Eqs.  (18)  and  (19)  that  Q'(xq)  “  0.  Thus,  Eq. 
(23)  reduces,  for  the  i  *  0  mode,  to 


Q(x0J 


L-  Qr  /2J1/2  eJK • 


(24) 


Equation  (23)  is  anologous  to  the  dispersion  relation  obtained  by  Ganguli 

and  Bakshi  (1982).  In  Eq.  (23)  the  nonlocal  effects  are  contained  in  the  two 

terms  on  the  right  hand  side,  one  of  which  contains  an  explicit  shear 

dependent  term  and  the  other  contains  a  shear  independent  nonlocal  term. 

This  equation  shows  that  in  the  p  /L  +  0  limit  the  term 

s  s 

containing  Q'  could  significantly  modify  the  growth  rate.  In  fact,  this 
shear  independent  nonlocal  term  drastically  reduces  the  growth  rate  in  the 
collisionless  kinetic  case.  Whereas,  Eq.  (24)  shows  that  in  the  fluid  case 
any  stabilizing  infuelnce  by  the  magnetic  shear  is  due  only  from  the  shear 
dependent  term  whiche  indicates  that  large  magnetic  shear  is  probably  needed 


to  stabilize  the  instability.  In  fact,  the  results  given  in  the  next 
section  support  this  observation.  Equation  (24)  is  solved  to  obtain  the 
eigenfrequencles  and  the  results  are  presented  in  the  next  section. 


III.  RESULTS 

We  solve  the  dispersion  equation  given  by  Eq.  (24)  numerically  in 
various  parameter  domains.  These  results  are  verified  by  directly  solving 
Eq.  (14)  using  a  numerical  shooting  code.  We  present  the  former  results  in 
the  following.  In  Fig.  1  we  plot  the  growth  rate  normalized  to  the  ion 
cyclotron  frequency  (y/0^)  versus  the  normalized  wavenumber  (k^Pg)  for 


Vs 


so,  v /a 

e  e 


0.01,  and  v^/Q^ 


0.01. 


treated  p  /L  as  an  external  parameter. 

8  S 


In  plotting  these  results  we 
The  solid  line  represents  t' 


shear-free  case  and  the  dashed  line  represents  the  case  where  p  /L  ■  0.01 

s  s 

We  see  from  the  figure  that  modes  with  small  k^p  s  are  not  strongly  affec 

by  the  magnetic  shear.  The  growth  rate  for  k^pg  “  0*3  reduces  from  0.24  ■«* 

0.20  for  p  /L  0.01,  whereas  the  growth  rate  for  k  p  »  0.7  reduces  by 
3  3  y  3 

about  30%  from  0.75  to  0.55. 


In  order  to  throw  some  light  on  the  magnetic  shear  required  to 

significantly  reduce  the  growth  rate,  we  plot  in  Fig.  2,  the  normalized 

growth  rate  versus  the  normalized  shear  length  (p  /L  ),  treating  p  /L  as 

s  s  s  s 

an  external  parameter.  We  choose  k  p  *0.3,  and  0.5  and  plot  the  growth 

y  s 

rate  for  the  parameters  of  Fig.  1.  This  figure  shows  that  moderate  to 

strong  magnetic  shears,  p  /L  >  0.01,  are  required  to  reduce  the  growth 

s  s 

rate  significantly.  For  example,  for  k  pg  ■  0.3  the  growth  rate  is  reduced 

by  about  50%  for  p  /L  ■  0.04  and  the  mode  Is  stabilized  for  p  /L  >  0.065, 

s  s  s  s 

and  for  k  p  »  0.5  the  growth  rate  is  reduced  by  50%  for  p  /L  =  0.04  but 
y  s  s  s 

the  mode  is  stabilized  for  P_/L_  -  0.078.  However,  the  theory  begins  to 


s  s 

break  down  for  large  shear  value  p  / L  >0.1. 

s  s 


We  find  that  for  large 


shears,  the  growth  rate  drops  linearly  as  a  function  of  inverse  shear 

-4 

length.  For  very  small  shears,  pg/Lg  *  10  the  growth  rate  remains 

essentially  constant  and  the  growth  rate  decreases  noticeably  for  p  / L  > 

s  s 

0.01.  In  an  ionospheric  environment  B  =  0.5  Gauss  and  pg  =*  2.5  meters,  and 

Eqs.  (1)  and  (3)  show  that  p  /L 

s  s 

0.1  A/m  and  1  A/m  ,  respectively. 


0.01  and  1.0  correspond  to  currents  of 


In  Fig.  3  we  present  the  wave  eigenfunctions  for  the  following 

-2  -2 

parameters:  V.  ■  50  c  ,  v  /a  ■  10  ,  v./Q.  -  10  and  for  k  p  ■  0.3, 

d  s  ei  e  ’  i  i  y  s 

Y/n  -  0.24.  The  solid  line  represents  the  real  part  of  the  wave  function 

($)  and  the  dashed  line  |$j.  Two  important  features  are  to  be  noted  in 

this  figure.  One  feature  is  that  the  wavepacket  localizes  around  xQ/pg  «  1> 

i.e.,  for  small  shear  lengths  k  /k  is  «  1,  indicating  that  the  mode  is 

z  y 

almost  perpendicularly  propagating.  This  can  also  be  seen  from  Eq.  (19) 
which  yields 

k  /k  -  2(io/k  V,J(p  /L  )  «  1  for  p  /L  -  10~7. 
z  y  y  d  v  s  s'  ss 

The  second  feature  is  that  the  width  of  the  wave  packet  is  of  the  order 

200  p  suggesting  a  localization  region  of  500  meters  for  p  =2.5  meters 

3  +  S 
(corresponding  to  NO  ions  in  the  ionosphere). 

The  effect  of  self-consistent  magnetic  shear  is  understood  by  plotting 

the  normalized  growth  rate  versus  the  normalized  drift  velocity.  For  this 

purpose,  we  use  Eq.  (4)  to  express  the  shear  length  in  terms  of  the  drift 

velocity  and  solve  Eq.  (24).  In  figure  4,  we  plot  the  growth  rate 

for  a  -  10  6,  10  5 ,  and  10  *  (uj2  /u>2  =*  102,  103,  and  104,  respectively; 

pe  ce 

c  »  500  m/s)  and  for  v  /ft  *  0.01,  v  /ft  *  0.01,  and  k  p  =  0.3.  We  also 
3  e  e  I  l  y  s 

give  the  shear-free  local  growth  rate  (curve  A)  to  compare  with  the  growth 
rate  with  shear.  Three  points  are  to  be  noted  in  this  figure:  (1)  a  *  10  ^ 
corresponds  to  ionospheric  parameters.  We  find  that  the  growth  rate  is  not 
too  different  from  the  shear-free  case.  In  fact,  the  growth  rate  curve 
overlaps  curve  A  and  eventually  turns  around  for  V^/cg  >  500.  (2)  Curve  B, 
which  gives  the  growth  rate  for  a  -  10  3,  shows  that  the  growth  rate  drops 

significantly  beyond  drift  velocities  >  150  c  ;  and  (3)  curve  C,  which 

-4 

represents  the  growth  rate  for  a  ^  10  ,  shows  that  the  optimum  drift 

velocity  (the  drift  velocity  for  which  the  growth  rate  is  a  maximum),  80  cg, 
is  much  smaller  for  this  case.  This  leads  to  the  conclusion  that  the 

optimum  drift  velocity  decreases  as  the  parameter  a  is  increased.  At  large 

drift  velocities,  the  self-consistent  magnetic  shear  produced  is  large 

enough  that  the  mode  is  stabilized  and  the  growth  rate  tends  to  zero. 
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IV.  *■  SCUSSION 


We  now  apply  the  results  to  the  high  latitude  Ionospheric  E  region.  We 

choose  v  /S)  «  10  v  /8  -  10  ^  and  V.  *•  20  c  .  For  V.  -  20  c_  the 

11  ee  ’  d  s  a  s 

corresponding  shear  length  Is  about  1000  km  (N0+  Ions,  n  -  IQ5  cm-3).  We 

find  that  the  growth  rate  Is  0.9  s-*,  for  modes  with  wavelengths  of  50 

meters  (k  p  *  0.3,  and  p  ■  2.5  meters),  which  is  comparable  with  the  local 
y  s  s 

growth  rate.  From  the  plots  of  the  wave  functions  we  see  that  the  mode  is 

localized  In  the  north-south  direction  In  a  region  of  Ax  **  200  p  *  500  m. 

s 

Moreover,  we  find  that  k  »  k  (k  -  1/Ax)  which  indicates  the  strong  two 

y  x  x 

dimensional  structure  of  the  mode  in  the  plane  perpendicular  to  the  magnetic 
field  during  the  linear  phase  of  the  instability. 

We  conclude  that  the  magnetic  shear  corresponding  to  normal  auroral 

2 

conditions  with  field  aligned  currents  of  the  order  of  pA/m  is  small 

(Lg  '  1000  km)  and  thus  does  not  apprecibly  alter  the  current  driven 

collisional  ion  cyclotron  mode  discussed  by  Chaturvedi  (1976).  Even  strong 

2 

currents  under  some  disturbed  conditions,  *»  135  pA/m  ,  reported  by  Burke  et 
al.  (1983)  do  not  produce  significant  shear  (Lg  *»  300  km)  in  the  magnetic 
field  to  have  any  effect  on  the  CIC.  Currents  of  the  order  mA/m^  or  more 
possibly  produce  significant  effects. 

Finally,  we  discuss  the  limitations  of  the  present  theory.  (1)  The 

mode  structure  equation  (Eq.  18)  and  the  subsequent  dispersion  relation  (Eq. 

23)  were  derived  under  the  assumption  that  kg/ ky  <  1  based  on  the  premise 

that  the  ion  parallel  motion  [V»(V^)  term  in  the  ion  continuity  equation] 

may  not  be  important.  Thus  the  results  presented  do  not  contain  ion 

parallel  motion  effects.  However,  we  extended  the  theory  to  include  the  ion 

parallel  motion  and  found  no  significant  changes  in  the  way  the  magnetic 

shear  affects  the  mode.  The  local  growth  rate  (in  the  limit  p  /L  **  0  ) 

s  s 

with  the  ion  parallel  motion  is  smaller  by  about  10 Z;  the  growth  rate  with 

magnetic  shear  is  consistently  smaller  but  has  similar  behaviour  as  shown 

earlier.  (2)  For  larger  shears,  p  /L  '0.1,  the  parabolic  expansion  of  the 

s  s 

potential  function  (Eq.  23),  used  to  derive  the  analytical  dispersion 

relation,  may  not  be  adequate.  This  is  due  to  the  fact  that  the  wave 

function  spreads  out  and  samples  non-parabolic  part  of  the  potential. 

Numerical  solutions  of  (Eq.  18)  do  confirm  the  spreading  of  the  wave 

packet.  Furthermore,  the  wave  packet  localizes  in  the  region  where  the 

effective  (k_/k„)  '1  because  magnetic  shear  introduces  an  effective  k_.  (3) 

z  y  z 
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We  have  done  an  MHD  analysis  and  arrived  at  the  mode  structure  equation  with 
a  potential  term,  Q(x)  [Eq.  (18)],  whose  derivative  vanishes  at  xQ,  whereas, 
in  the  collisionless  case  the  derivative  of  Q(x)  at  xQ  is  finite,  thus 
drastically  affecting  the  collisionless  kinetic  ion  cyclotron  instability 
(Ganguli  and  Bakshi,  1982).  To  further  examine  the  above  three  aspects  we 
will  present  the  kinetic  analysis  of  the  collisonal  ion  cyclotron  waves  in 
the  presence  of  magnetic  shear  in  a  future  paper.  For  higher  drift 
velocities  (comparable  to  the  E  *  B  drift  velocities)  where  the  self- 
consistent  magnetic  shear  plays  a  significant  role  in  determining  the  mode 
structure,  a  complete  treatment  is  needed  which  includes  an  electric  field 
along  the  y-direction  and  the  related  velocity  shear. 

In  conclusion,  we  have  examined  the  influence  of  magnetic  shear  on  the 

collisional  current  driven  ion  cyclotron  instability.  Self  consistent 

magnetic  shear  corresponding  to  moderate  drift  velocities  (V^  ■'0.5  -  5 

km/s)  near  the  threshold  velocity  for  the  instability  does  not  have 

significant  effects  on  the  instability.  We  find  that  the  mode  is  almost 

perpendicularly  propagating  and  is  localized  in  a  region  extending  upto  a 

few  hundred  kilometers  in  the  direction  of  the  magnetic  field  under  auroral 

conditions.  However,  in  domains  where  the  plama  density  is  high  such  that 
2  2 

the  parameter  w  /uce  Is  large,  we  find  that  the  growth  rate  of  a 
particular  mode  maximizes  at  an  optimum  drift  velocity  much  larger  than  the 
sound  speed.  Finally,  we  find  that  strong  shears  (ps/Lg  >  0.05)  signifi¬ 
cantly  reduce  the  growth  rate  and  stabilize  the  collisional  current  driven 
ion  cyclotron  instability.  This  strong  shear  corresponds  to  large  parallel 
drift  velocities  (V^  »  cg) ,  as  seen  in  fig.  4. 
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Figure  1.  Plot  of  normalized  growth  rate  (y/0^)  versus  the  normalized 

wave  number  (k  p  ).  The  parameters  used  are  v  /O  *  0.01, 
y  8  6  6 

Vf/Oi  ■  0.01,  and  V^/Cg**  50.  The  solid  line  represents  the 
shear  free  case,  and  the  dashed  line  represents  the  case  where 
p  /L  ■  0.01. 

s  s 
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Figure  2.  Normalized  growth  rate  (y/Q.)  versus  p  /L  treating  p  /L  as  an 

X  3  3  3  8 

external  parameter  for  V^/Cg  ■  50,  ly>a  ■  0.3,  and  0.5,  and  for 

v  /0  ■  0.01,  and  v./O,  ■  0.01. 

e  e  11 
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Figure  4.  Growth  rate  versus  the  drift  velocity  curves  with  and  without 
self-consistent  magnetic  shear  for  the  following  parameters: 

Curve  A  represents 


k  p 
y  s 


0.3,  V  /n 
e  e 


0.01,  and  v./o,  ■  0.01. 

—6*  * 

the  shear-free  case  and  a  m  10  .  Curves  B  and  C  are  for  a  « 

e  2  2 

10  and  10  4,  respectively,  and  these  correspond  to  m _ /<•) 

.  1  pe  ce 

10^  and  10  . 
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HANSCOM  AFB,  MA  01731 
OICY  ATTN  J.  DEAS 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DI VISION/YSEA 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01732 
OICY  ATTN  YSEA 


HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION/DC 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731 

01CY  ATTN  DCKC  MAJ  J.C.  CLARK 

COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRIGHT-PATTERSON  AFB,  OH  45433 
01CY  ATTN  NICD  LIBRARY 

01CY  ATTN  ETDP  B.  BALLARD 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIFFISS  AFB,  NY  13441 

01CY  ATTN  DOC  LIBRARY/TSLD 
01CY  ATTN  OCSE  V.  COYNE 

SAMSO/SZ 

POST  OFFICE  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
(SPACE  DEFENSE  SYSTEMS) 

01CY  ATTN  SZJ 

STRATEGIC  AIR  COMMAND/XPFS 
OFFUTT  AFB,  NB  68113 

01CY  ATTN  ADWATE  MAJ  BRUCE  BAUER 
01CY  ATTN  NRT 

01CY  ATTN  DOK  CHIEF  SCIENTIST 

SAMSO/SK 
P.O.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

01CY  ATTN  SKA  (SPACE  COMM  SYSTEMS) 
M.  CLAVIN 

SAMSO/MN 

NORTON  AFB,  CA  92409 
(MINUTEMAN) 

01CY  ATTN  MNNL 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01731 

01CY  ATTN  EEP  A.  LORENTZEN 

DEPARTMENT  OF  ENERGY 
LIBRARY  ROOM  G-042 
WASHINGTON,  D.C.  20545 

01CY  ATTN  DOC  CON  FOR  A.  LABOWITZ 


DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.O.  BOX  5400 
ALBUQUERQUE,  NM  87115 

01CY  ATTN  DOC  CON  FOR  D.  SHERWOOD 

EG&G,  INC. 

LOS  ALAMOS  DIVISION 

P.O.  BOX  809 

LOS  ALAMOS,  NM  85544 

01CY  ATTN  DOC  CON  FOR  J.  BREEDLOVE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 

01CY  ATTN  DOC  CON  FOR  TECH  INFO  DEPT 

01CY  ATTN  DOC  CON  FOR  L-389  R.  OTT 

01CY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 

01CY  ATTN  DOC  CON  FOR  L-46  F.  SEWARD 

LOS  ALAMOS  NATIONAL  LABORATORY 

P.O.  BOX  1663 

LOS  ALAMOS,  NM  87545 

01CY  ATTN  DOC  CON  FOR  J.  WOLCOTT 

01CY  ATTN  DOC  CON  FOR  R.F .  TASCHEK 

01CY  ATTN  DOC  CON  FOR  E.  JONES 

01CY  ATTN  DOC  CON  FOR  J.  MALIK 

01CY  ATTN  DOC  CON  FOR  R.  JEFFRIES 

01CY  ATTN  DOC  CON  FOR  J.  ZINN 

01CY  ATTN  DOC  CON  FOR  P.  KEATON 

01CY  ATTN  DOC  CON  FOR  D.  WESTERVELT 
01CY  ATTN  D.  SAPPENFIELD 

SAND1A  LABORATORIES 
P.O.  BOX  5800 
ALBUQUERQUE,  NM  87115 

01CY  ATTN  DOC  CON  FOR  W.  BROWN 

OICY  ATTN  DOC  CON  FOR  A.  THORNBROUGH 

OiCY  ATTN  DOC  CON  FOR  T.  WRIGHT 

OICY  ATTN  DOC  CON  FOR  D.  DAHLGREN 

OICY  ATTN  DOC  CON  FOR  3141 

OICY  ATTN  DOC  CON  FOR  SPACE  PROJECT  DIV 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  B.  MURPHEY 

OICY  ATTN  DOC  CON  FOR  T.  COOK 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  CON  DR.  YO  SONG 
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OTHER  GOVERNMENT 

DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  D.C.  20234 

OICY  (ALL  CORRES:  ATTN  SEC  OFFICER  FOR) 

INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  &  INFO  ADMIN 
BOULDER,  CO  80303 

01CY  ATTN  A.  JEAN  (UNCLASS  ONLY) 

01CY  ATTN  W.  UTLAl'T 
01CY  ATTN  D.  CROMBIE 
01CY  ATTN  L.  BERRY 

NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 
01CY  ATTN  R.  GRUBB 
01CY  ATTN  AERONOMY  LAB  G.  REID 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
01CY  ATTN  I.  GARFUNKEL 
01CY  ATTN  T.  SALMI 
01CY  ATTN  V.  JOSEPHSON 
01CY  ATTN  S.  BOWER 
01CY  ATTN  D.  OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

01CY  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOC.,  INC. 

1901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 
01CY  ATTN  L.  SLOAN 
OICY  ATTN  R.  THOMPSON 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY,  CA  94701 
OICY  ATTN  J.  WORKMAN 
OICY  ATTN  C.  PRETTIE 
OICY  ATTN  S.  BRECHT 


BOEING  COMPANY,  THE 
P.O.  BOX  3707 
SEATTLE,  WA  98124 
OICY  ATTN  G.  KEISTER 
OICY  ATTN  D.  MURRAY 
OICY  ATTN  G.  HALL 
OICY  ATTN  J.  KENNEY 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 
OICY  ATTN  D.B.  COX 
OICY  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  20734 
OICY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

OICY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
OICY  ATTN  H.  LOGSTON 
OICY  ATTN  SECURITY  (PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  Ullshlre  Blvd. 

Santa  Monica,  Calif  90401 
OICY  ATTN  C.S.  GABBARD 

ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
OICY  ATTN  J.  ROBERTS 
OICY  ATTN  JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 

OICY  ATTN  M.H.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NY  13201 
OICY  ATTN  F.  REIBERT 


GENERAL  ELECTRIC  TECH  SERVICES  CO.  INC. 
RMES 

COl'RT  STREET 
SYRACUSE,  NY  13201 
01CY  ATTN  G.  MILLMAN 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
OICY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
01CY  ATTN  TECHNICAL  LIBRARY 
OICY  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

CTE  SYL VANIA,  INC. 

ELECTRONICS  SYSTEMS  CRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

OICY  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

OICY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

OICY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  VA  22311 
OICY  ATTN  J.M.  AEIN 
OICY  ATTN  ERNEST  BAUER 
OICY  ATTN  HANS  WOLFARD 
OICY  ATTN  JOEL  BENGSTON 

INTL  TEL  6  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 
OICY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

11011  TORREYANA  ROAD 
P.O.  BOX  85154 
SAN  DIEGO,  CA  92138 
OICY  ATTN  J.L.  SPERLING 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

OICY  ATTN  DOCUMENT  LIBRARIAN 
OICY  ATTN  THOMAS  POTEMRA 
OICY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
OICY  ATTN-  T.  MEAGHER 

KAMAN  TEMPO  -CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
OICY  ATTN  DASIAC 
OICY  ATTN  WARREN  S.  KNAPP 
OICY  ATTN  WILLIAM  MCNAMARA 
OICY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
OICY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC. 

3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-12 
OICY  ATTN  W.L.  IMHOF  DEPT  52-12 
OICY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
OICY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 
OICY  ATTN  R.  HEFFNER 


M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  HA  02173 

OICY  ATTN  DAVID  M.  TOWLE 
OICY  ATTN  L.  LOUGHLIN 
OICY  ATTN  D.  CLARK 
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MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01CY  ATTN  N.  HARRIS 
01CY  ATTN  J.  MOULE 
01CY  ATTN  GEORGE  MROZ 
01CY  ATTN  W,  OLSON 
01CY  ATTN  R.W.  HALPRIK 
01CY  ATTN  TECHNICAL  LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 

735  STATE  STREET 

SANTA  BARBARA,  CA  93101 


01CY 

ATTN 

P.  FISCHER 

01CY 

ATTN 

W.F.  CREVIER 

01CY 

ATTN 

STEVEN  L.  GUTSCHE 

01CY 

ATTN 

R.  BOGUSCH 

01CY 

ATTN 

R.  HENDRICK 

01CY 

ATTN 

RALPH  KILB 

01CY 

ATTN 

DAVE  SOWLE 

01CY 

ATTN 

F.  FA JEN 

01CY 

ATTN 

M.  SCHEIBE 

01CY 

ATTN 

CONRAD  L.  LONGMIRE 

01CY 

ATTN 

B.  WHITE 

MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.E. 
ALBUQUERQUE,  NEW  MEXICO  87106 
01CY  R.  STELLINGWERF 
01CY  M.  ALME 
01CY  L.  WRIGHT 

MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

01CY  ATTN  JOHN  MORGANSTERN 
01CY  ATTN  G.  HARDING 
01CY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
01CY  ATTN  W.  HALL 
01CY  ATTN  W.  FOSTER 


PENNSYLVANIA  STATE  UNIVERSITY 

IONOSPHERE  RESEARCH  LAB 

318  ELECTRICAL  ENGINEERING  EAST 

UNIVERSITY  PARK,  PA  16802  *  < 

(NO  CLASS  TO  THIS  ADDRESS) 

01CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

01CY  ATTN  IRVING  L.  KOFSKY  I  ( 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  UA  98009 

01CY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC.  . 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  &  D  ASSOCIATES 

P.O.  BOX  9695 

MARINA  DEL  REY,  CA  90291 


01CY 

ATTN 

FORREST  GILMORE 

1 

( 

01CY 

ATTN 

WILLIAM  B.  WRIGHT,  JR. 

01CY 

ATTN 

ROBERT  F.  LELEVIER 

01CY 

ATTN 

WILLIAM  J.  KARZAS 

01CY 

ATTN 

H.  ORY 

01CY 

ATTN 

C.  MACDONALD 

01CY 

ATTN 

R.  TURCO 

01CY 

ATTN 

L.  DeRAND 

* 

( 

01CY 

ATTN 

W.  TSAI 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
01CY  ATTN  CULLEN  CRAIN 
01CY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

01CY  ATTN  BARBARA  ADAMS 


PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

01CY  ATTN  E.C.  FIELD,  JR. 


RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

01CY  ATTN  VINCE  TRAPANI 


i 
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SCIENCE  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

OICY  ATTN  LEVIS  M.  LINSON 
OICY  ATTN  DANIEL  A.  HAMLIN 
01CY  ATTN  E.  FRIEMAN 
01CY  ATTN  E.A.  STRAKER 
OICY  ATTN  CURTIS  A.  SMITH 
01CY  ATTN  JACK  MCDOUGALL 

SCIENCE  APPLICATIONS,  INC 
1710  GOODRIDGE  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 

OICY  ATTN  DONALD  NEILSON 
OICY  ATTN  ALAN  BURNS 
OICY  ATTN  G.  SMITH 
OICY  ATTN  R.  TSUNODA 
OICY  ATTN  DAVID  A.  JOHNSON 
OICY  ATTN  WALTER  G.  CHESNUT 
OICY  ATTN  CHARLES  L.  RINO 
OICY  ATTN  WALTER  JAYE 
OICY  ATTN  J.  VICKREY 
OICY  ATTN  RAY  L.  LEADA8RAND 
OICY  ATTN  G.  CARPENTER 
OICY  ATTN  G.  PRICE 
OICY  ATTN  R.  LIVINGSTON 
OICY  ATTN  V.  GONZALES 
OICY  ATTN  D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 
OICY  ATTN  W.P.  BOQUIST 

TOYON  RESEARCH  CO. 

P.O.  Box  6890 
SANTA  BARBARA,  CA  93111 
OICY  ATTN  JOHN  ISE,  JR. 

OICY  ATTN  JOEL  GARBARINO 


TRW  DEFENSE  &  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
OICY  ATTN  R.  K.  PLEBUCH 
OICY  ATTN  S.  ALTSCHULER 
OICY  ATTN  D.  DEE 
OICY  ATTN  D/  STOCKWELL 
SNTF/1575 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
OICY  ATTN  W.  REIDY 
OICY  ATTN  J.  CARPENTER 
OICY  ATTN  C.  HUMPHREY 


IONOSPHERIC  MODELING  DISTRIEVTION  LIST 
(UNCLASSIFIED  0:."Li) 

PLEASE  DISTRIBUTE  ONE  COPY  TO  EACH  OF  THE  FOLLOWING  PEOPLE  (UNLESS  OTHERWISE 
NOTED) 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 

Dr.  P.  MANGE  -  CODE  A 101 
Dr.  P.  GOODMAN  -  CODE  4180 

A.F.  GEOPHYSICS  LABORATORY 
L.G.  HANSCOM  FIELD 
BEDFORD,  MA  01730 
DR.  T.  ELKINS 
DR.  W.  S WIDER 
MRS.  R.  SAGALYN 
DR.  J.M.  FORBES 
DR.  T.J.  KENESHEA 
DR.  W.  BURKE 
DR.  H.  CARLSON 
DR.  J.  JASPERSE 

BOSTON  UNIVERSITY 
DEPARTMENT  OF  ASTRONOMY 
BOSTON,  MA  02215 
DR.  J.  AARONS 

CORNELL  UNIVERSITY 
ITHACA,  NY  1A850 
DR.  W.E.  SWARTZ 
DR.  D.  FARLEY 
DR.  M.  KELLEY 

HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MA  02138 
DR.  M.B.  HcELROY 
DR.  R.  LINDZEN 


NASA 

GODDARD  SPACE  FLICHT  CENTER 
GREENBELT ,  MD  20771 
DR.  K.  MAEDA 
DR.  S.  CURTIS 
DR.  M.  DUBIN 

DR.  N.  MAYNARD  -  CODE  696 
COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 
DR.  T.  CZUBA 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  CA  92152 

MR.  R.  ROSE  -  CODE  5321 

NOAA 

DIRECTOR  OF  SPACE  AND 
ENVIRONMENTAL  LABORATORY 
BOULDER,  CO  80302 
DR.  A.  GLENN  JEAN 
DR.  G.U.  ADAMS 
DR.  D.N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.F.  DONNELLY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VA  22217 
DR.  G.  JOINER 


INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY/NAVY  DRIVE 
ARLINGTON,  VA  22202 
DR.  E.  BAUER 

MASSACHUSETTS  INSTITUTE  OF 
TECHNOLOGY 

PLASMA  FUSION  CENTER 
LIBRARY,  NW16-262 
CAMBRIDGE,  MA  02139 


PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PA  16802 
DR.  J.S.  NISBET 
DR.  P.R.  ROHRBAUGH 
DR.  L.A.  CARPENTER 
DR.  M.  LEE 
DR.  R.  DIVANY 
DR.  P.  BENNETT 
DR.  F.  KLEVANS 

SCIENCE  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 
DR.  D.A.  HAMLIN 
DR.  E.  FRIEMAN 


STANFORD  UNIVERSITY 
STANFORD,  CA  94305 
DR.  ?.M.  BANKS 

U.S.  ARMY  ABERDEEN  RESEARCH 
AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  MD 

DR.  J.  HEIMERL 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 
DR.  L.E.  LEE 

UNIVERSITY  OF  CALIFORNIA, 
BERKELEY 

BERKELEY,  CA  94720 
DR.  M.  HUDSON 

UNIVERSITY  OF  CALIFORNIA 
LOS  ALAMOS  SCIENTIFIC  LABORATORY 
J-10,  MS-664 
LOS  ALAMOS,  NM  87545 
DR.  M.  PONGRATZ 
DR.  D.  SIMONS 
DR.  G.  BARASCH 
DR.  L.  DUNCAN 
DR.  P.  BERNHARDT 
DR.  S.P.  GARY 

UNIVERSITY  OF  MARYLAND 
COLLECE  PARK,  MD  20740 
DR.  K.  PAPADOPOULOS 
DR.  E.  OTT 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 
DR.  R.  GREENWALD 
DR.  C.  MENG 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 
DR.  N.  ZABUSKY 
DR.  M.  BIONDI 
DR.  E.  OVERMAN 


UNIVERSITY  OF  TEXAS 
AT  DALLAS 

CENTER  FOR  RESEARCH  SCIENCES 
P.O.  BOX  688 
RICHARDSON,  TX  75080 
DR.  R.  HEELIS 
DR.  W.  HANSON 
DR.  J.P.  McCLURE 

UTAH  STATE  UNIVERSITY 
4TH  AND  8TH  STREETS 
LOGAN,  UTAH  84322 
DR.  R.  HARRIS 
DR.  K.  BAKER 
DR.  R.  SCHUNK 
DR.  J.  ST. -MAURICE 

PHYSICAL  RESEARCH  LABORATORY 
PLASMA  PHYSICS  PROGRAMME 
AHMEDABAD  380  009 
INDIA 

P.J.  PATHAK,  LIBRARIAN 

LABORATORY  FOR  PLASMA  AND 
FUSION  ENERGY  STUDIES 
UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20742 
JHAN  VARY AN  HELLMAN, 
REFERENCE  LIBRARIAN 
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